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CHROMATOGRAFHIC  THEORY 


Gas  chromatography  is  a  versatile,  relatively  new  analytical  tech- 
nique.    Chromatography's  history  dates  back  to  the  work  of  F.  F.  Runge  in 
1850  (1).     The  usefulness  of  gas  chromatography  became  apparent  in  the 
decade  beginning  in  1950.    Since  then  researchers  have  found  gas  chroma- 
tography useful  in  a  wide  range  of  separation  and  identification  problems 
involving  both  organic  and  inorganic  problems. 

Gas -liquid  chromatography  is  a  process  based  on  partition  equilib- 
rium of  a  vapor  between  a  stationary  liquid  phase  dispersed  on  a  porous 
solid  support  and  a  moving  gas  phase.    Volatile  mixtures  can  be  separated 
into  components  when  each  component  has  a  different  distribution  coeffi- 
cient between  the  liquid  phase  and  gas  phase . 

The  time  between  injection  of  a  compound  into  the  column  and  its 
elution  and  detection  is  t^,  the  experimental  retention  time.    This  reten- 
tion time  is  directly  proportional  to  the  volume  of  gas  that  is  needed  to 

elute  the  compound.    This  volume  is  referred  to  as  V^,  the  experimental 

R 

retention  voliame.    These  parameters  are  related  by 

V  =  \  '  (1) 

where  v  is  the  carrier  gas  flow  rate  measured  at  the  end  of  the  column. 

The  ratio  of  the  amount  of  solute  in  the  liquid  phase  to  that  in  the 
gas  phase  is  denoted  by  k,  commonly  called  the  capacity  ratio.    In  gas- 
liquid  chromatography,  the  retention  volume,  the  partition  ratio,  and  the 
retention  volume  of  some  nonsorbed  gas  V  ,  are  related  by 

VR=V^(l+k).  (2) 

We  may  eliminate  the  dead  volume,  V  ,  of  the  column 
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which  is  ^  constant  contribution  to  retention  volume  by  defining  a  net 
retention  volume  V ' 


R" 


Similarly  we  may  define  the  net  retention  time,  t'  : 

R 

*'R=V^a  =  ^\- 

in  which  t^  is  the  time  required  for  elution  of  a  nonsorbed  compound. 

The  flow  rate  of  the  carrier  gas  measured  at  the  end  of  the  column 

can  be  related  to  the  velocity  of  the  carrier  gas  through  the  column-. 

However,  the  linear  velocity  is  not  constant;  it  varies  as  a  function  of 

position  along  the  column.    This  variation  is  a  function  of  the  pressure 

drop  which  can  be  measured  at  points  along  the  column.    The  pressure  at 

position  X  is  related  to  the  inlet  and  outlet  pressure  by 

X      (P,/P^)^-(P,/P,)^  . 

-=  1   (5) 

L      (Pi/Po)'  -1 

where  L  is  the  length  of  the  column,  p    and  p    are  inlet  and  outlet  pres- 

i  o 

wures  respectively,  and  p^  is  the  pressure  at  position  x,  where  x  is  the 
distance  measured  from  the  column  inlet.    The  above  equation  is  referred 
to  as  the  pressure  correction  factor  and  was  developed  by  James  and 
Martin. 

The  linear  velocity  of  gas  at  position  x  may  be  calculated  by  sub- 
stitution of  the  relation 

V^x=Px/Po  (6) 

into  equation  number  (5),  where  u^  is  the  velocity  at  the  outlet  of  the 
column  and  u^  is  the  velocity  at  position  k.    It  is  useful  to  know  the 
average  pressure,  p: 
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Pdx 

1—  '  (7) 


This  integration  gives 

-  .  p  =  2  p  (Pi/Po)^-i .  ' 

^       (Pi/Po)  -1 

The  average  velocity  u  along  the  column  may  he  expressed  as  (2) : 

^  -  -I—  •  (9) 

P 

The  capacity  ratio  is  related  to  the  partition  coefficient,  K;  by 

k  =  K  V  /V  (10) 
1    g  ' 

where  K  is  the  concentration  in  the  liquid  phase  divided  by  the  concen- 
tration in  the  gas  phase.    The  letters  V    and  V    represent  the  volumes  of 
liquid  and  gas  phases  contained  in  unit  length  of  the  column.    The  par- 
tition coefficient,  K,  is  a  thermodynamic  parameter  which  can  be  used  to 
identify  a  sample  component,  if  the  partition  coefficient  is  directly 
proportional  to  the  net  retention  time,  t'^.     This  can  be  shown  by  solving 
equation  {h)  for  the  partition  coefficient  and  expressing  t    as  l/u: 


=  L/-[K(V3_/VJ].  _  (11) 


The  partition  of  a  component  in  the  column  between  the  liquid  phase 
on  the  inert  support  and  the  mobile  gas  phase  may  be  treated  as  a  chemical 
equilibrium.     The  equilibrium  constant  is  the  partition  coefficient;  which 
is  the  ratio  of  the  concentration  of  solute  in  the  liquid  phase  to  the 


concentration  of  solute  in  the         phase.    The  partition  coefficient  is 
a  function  of  temperatiire  and  of  the  molar  free  energy  charge,  AF.  Thus 
the  distribution  coefficient,  K,  can  also  be  related  to  the  molar  heat  of 
solution  and  molar  entropy: 

InK  =  AF/rT  =  -Ah/RT+-AS/r  (12) 

where  T  is  absolute  temperature,  and  R  is  the  molar  gas  constant,  AH  is 
the  molar  heat  of  solution  and  AS  is  the  molar  entropy.    From  equation 
(12)  it  is  clear  that  K,  and  any  retention  parameter  based  on  K,  depend 
strongly  on  the  temperature.    It  can  be  observed  from  equation  (12)  that 
an  increase  in  temperature  will  decrease  the  value  of  the  distribution 
coefficient  and  thus  decrease  the  retention  time. 

Operation  at  high  column  temperature  will  shorten  analysis  time, 
but  this  is  achieved  at  the  expense  of  resolution  for  low-boiling  com- 
ponents.   "When  the  column  temperature  is  raised  during  a  run,  low-boiling 
components  can  be  eluted  with  adequate  resolution  at  the  low  initial  tem- 
perature, and  heavy  components  can  still  be  eluted  in  reasonable  time  at 
the  elevated  temperatures  reached  toward  the  end  of  the  run. 

In  addition  to  K,  other  parameters  such  as  mass  flow  rate  and 
average  velocity  are  affected  during  temperature -time  programs.  The 
change  in  temperature  affects  the  viscosity  of  the  carrier  gas,  which  in 
turn  affects  the  average  velocity  of  the  carrier.    At  constant  inlet  and 
outlet  pressures  the  velocity  profile  of  the  column  remains  constant;  that 
is,  at  any  point  on  the  coliimn  the  ratio  of  the  average  velocity  to  local 
velocity  is  independent  of  T.    Therefore,  only  the  average  velocity  need 
be  known  ais  a  function  of  temperature.       '  '   •  " 
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The  viscosity- temperature  relation  was  expressed  by  Golay  (3)  as: 

=  (TyT)Y  (13) 

vhere  t]^  is  the  viscosity  at  starting  temperature       temperature  in 
degrees  Kelvin^  and  "n^  is  the  viscosity  at  absolute  temperature  T.  For 
helium  y  =  0.^.    At  constant  flow  at  the  outlet  of  the  column,  the  average 
velocity  depends  upon  temperature,  and  its  variation  can  be  predicted 
from  the  gas  laws. 

At  constant  inlet  and  outlet  pressures  the  pressure  profile  along 
the  column  is  independent  of  temperature  only  if  the  column  is  uniformly 
heated. 

The  gas  flow  in  a  column  under  constant  inlet  and  outlet  pressure  is 
not  a  constant  during  a  programmed-temperature  run.    Therefore,  a  detector 
which  is  insensitive  to  gas  flow  is  req.uired.    Hydrogen  flame  ionization 
detectors  are  useful  for  such  systems,  because  they  are  largely  insensitive 
to  flow  rate .     •  .  ■ 

Although  the  distribution  coefficient  can  be  calculated  from  reten- 
tion data,  this  requires  knowledge  of  the  value  of  a  number  of  experimental 
parameters,  some  of  which  are  not  easily  accessible.     The  difficulties  may 
be  eliminated  in  isothermal  chromatography  by  the  use  of  relative  reten- 
tion data,  whereby  the  retention  is  expressed  relative  to  the  retention  of 
one  or  more  of  a  series  of  reference  compounds.     This  concept  has  been 
developed  by  Kovats  {k) .  '  '  " 

In  programmed  temperature  gas  chromatography,  retention  .data  may  also 
be  used  for  identification  purposes,  but  the  situation  is  complicated 
because  the  retention  of  a  given  compound  depends  on  the  time -temperature" 
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program  of  the  column  during  the  analysis.    Excellent  theoretical  work  has 
been  done  on  the  calculation  and  interpretation  of  retention  data  in  pro- 
grammed temperature  analysis,  notably  by  Golay  (5)  and  Giddings  (6),  but 
practical  applications  are  hampered  by  uncertainty  about  actual  column 
temperatures. 

A  retention  index  for  programmed  temperature  chromatography  which  is 
based  upon  elution  temperature  of  the  compound  has  been  proposed  by 
Guiochon  (7)-    The  nimerical  value  of  Guiochon's  index  is  approximately 
equal  to  the  Kovats  retention  index. 

Programmed  temperature  retention  data  can  have  a  large  error  due  to 
the  lack  of  reproducibility  in  the  temperature  prograj3i(8) .    Habgood  and 
Harris  greatly  increased  the  usefulness  of  the  retention  indices,  Ip  (x), 
in  programmed  temperature  gas  chromatography  by  showing  that  retention 
indices  can  be  calculated  from  isothermal  data  (9).    Thus  the  retention 
time  or  temperature  in  temperature  programmed  gas  chromatography  can  be 
predicted  from  Isothermal  data. 

The  usefulness  of  programmed  temperature  gas  chromatography  could  be 
significantly  increased  if  reproducible  and  accurate  linear  temperature 
programs  could  be  developed.    The  retention  indices  for  compounds  eluted 
from  a  column  with  a  known  linear  temperature  program  would  be  as  useful 
as  isothermal  retention  indices.     To  date,  there  are  no  commercially 
available  chromato graphs  which  have  temperature  programs  which  are  both 
accurate  and  precise  to  the  degree  needed  for  useful  programmed  tempera- 
ture retention  indices. 
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EESIGN  OF  THEEMOSTAT: 


I. 


Theory 


Presently  available  conimercial  equipment  depends  upon  heat  transfer 
from  a  thermostat  to  the  column  by  air  convection  or  direct  conduction. 
In  such  systems  there  is  an  inherent  time  lag  between  the  air  temperature 
and  the  column  temperature.    This  leads  to  an  uncertainty  in  the  column 
temperature,  the  magnitude  of  which  depends  on  the  time  lag.    Time  lag  'f 
can  be  calculated  from 


where  G  is  heat  transfer  coefficient  and  C  is  column  heat  capacity.  The 
time  lag  can.be  related  to  temperature  lag,  tsK ,  between  temperature  of 
the  thermostat  and  the  actual  temperature  of  the  column  by 


in  which  dT/dt  -is  the  rate  of  heating  of  the' thermostat.    The  minimtim  T' 
for  commercial  equipment  is  about  15  seconds.    It  may  be  as  large  as  100 
seconds  or  more,  but  then  the  equipment  is  seldom  claimed  as  suitable  for 
programmed  temperature  gas  chromatography.    Thus  the  magnitude  of  the 
temperature  difference  depends  upon  the  heating  rate;  for  a  normal  pro- 
gram rate,  the  temperature  difference  can  be  as  much  as  15° C  and  is  sel- 
dom less  than  1°C.     This  temperatxure  discrepancy  in  turn  will  affect  the 
reproducibility  of  the  program  and  thus  the  reproducibility  of  the  reten- 
tion data. 

The  time -temperature  discrepancies  can  be  corrected  by  placing  a 
thermocouple  on  the  column  and  measuring  the  time -temperature  program-  ^ 
for  each  run.    This  would  increase  the  usefulness  of  programmed 


^/^eff  =  " 


M  =  (dT/dt  jl' 


(15) 
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temperature  retention  indices,  but  also  complicate  a  now  simple  analytical 
method.    The  programmed  runs  would  be  simple  and  provide  useful  retention 
data  if  the  column  time -temperature  relation  were  accurately  known  and 
reproducible .    A  thermostat  has  been  developed  which  has  temperature  con- 
trol to  0.2°C  during  a  repeatable  program  run.    This  is  achieved  by  ohmic 
heating  of  an  adiabatically  mounted  column. 

The  purpose  of  this  thermostat  is  to  achieve  a  temperature  program 
of  a  gas  chromatographic  column  which  can  be  accurately  predicted  from 
constants  obtained  from  the  column.    Heat  is  generated  in  the  column  by 
means  of  an  electrical  current.    The  temperature  of  the  surrounding  air 
bath  is  kept  at  the  same  temperature  as  the  column.    Thus  there  is  no 
heat  transfer  to  or  from  the  column.     In  the  absence  of  heat  loss  the 
temperature  of  the  column  is  determined  from  the  electrical  energy  input 
and  the  heat  capacity  of  the  column.    Stray  heat  transferred,    A      to  the 
column  depends  on  two  factors:    the  heat  transfer  coefficient  G.  between 
the  column  and  its  surroundings,  and  the  temperature  gradient  AT  between 
the  column  and  its  surroundings ;  i . e . , 

^  "  ^TG^  ^  ,  ...   ^    .  (16) 

where  A  is  the  heat  leak  from  the  thermostat.    For  proper  operation  G^ 
and  AT  should  be  small.    We  can  recognize  two  factors  contributing  to 
G-j_:    air  convection,  and  conduction  through  column  connectors  and  tem- 
perature sensors.     Column  connectors  are  thermally  and  electrically 
insulated  by  small  Teflon  inserts  (Fig.  l) ;  thin  wire  thermocouples. 
(#36  copper- cons tantan)  are  used  as  temperature  sensors.    The  thermal 
gradient  AT  across  all  heat  leaks  is  kept  small  by  control  of  the  tem- 
perature of  the  air  bath  surrounding  the  column  container  and  by 
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providing  good  thermal  contact  between  the  air  bath  and  the  leads  and 
tubing  connected  to  the  column.    The  result  of  these  measures  can  be 
described  in  terms  of  an  effective  heat  transfer  coefficient  G^ff  from 
the  column  to  atmosphere . 

The  power  output,  P,  to  the  heater  of  the  thermostat  at  constant 
load  is  directly  proportional  to  the  temperature  difference  between  the 
column  and  the  thermostat. 

P  =  A'(T^-T^)+Po  ■  (IT) 

where  P  is  the  power,  A'  is  a  transfer  function  indicating  the  sensi- 
tivity of  the  control  system,  T    is  the  thermostat  temperature  and  P  is 
the  steady  input  power  to  the  thermostat.    The  temperature  difference 
between  thermostat  and  the  laboratory  is  proportional  to  P,  so  that  we 
may  write : 

\  V^o  "  ^^^c'^T^^Y^o     .      -     •    '  ■  ^^^^ 

where  y  is  a  power-to-temperature  conversion  factor,  and  T    is  the  labo- 

o 

ratory  air  temperature.    If  a  point  temperature  T.  is  defined  as 

^s  =  VY^o  '  '        ,  ■      "  (19) 

then  equation  (l8)  can  be  written:  .      .  ■ 

A  =  (Trj-Tg)/  (T^-Trp).  (20) 

The  heat  leak  away  from  the  colimn  to  the  thermostat  is  proportional 
to  the  heat  transfer  coefficient  and  the  temperature  gradient:      •  ' 

■     A  =  (T^-T^)"G. .    *  .  •  (21)- 
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As  seen  by  the  column,  the  heat  leak  may  also  be  written 


(22) 


Combining  ec[uation  (21)  and  (22)  and  simplifying  gives 


(23) 


If  equation  (22)  is  rearranged  and  T  +YP    substituted  for  T  ,  an 

o      o  s 

equation  for  a  line  of  the  form  Y  =  ax+b  is  obtained: 


where  AT  is  the  temperature  difference  between  the  column  and  the  labo- 
ratory, and  T    is  • 
P  o 

To  obtain  a  linear  time -temperature  program,  allowance  must  be  made 
for  changes  in  heat  capacity  and  electrical  resistance  of  the  column  with 
temperature.     The  resistance  of  the  column  increases  as  the  temperature 
increases  and  the  heat  capacity  decreases  as  the  temperature  increases. 
With  constant  voltage  input  to  the  column  the  two  act  in  opposite 
directions. 

A  differential  thermocouple  monitors  the  temperature  difference 
between  the  thermostat  and  the  column:  this  thermocouple  controls  the 
thermostat  heating  power  in  such  a  way  that  the  temperature  difference 
is  minimized  at  all  times.    The  control  loop  which  regulates  the  power  to 
the  thermostat  heaters  consists  of  a  Philbrick  voltage    preamplifier  and  a 
Honeywell  SCR  proportional  temperature  controller.    Details  are  given  on 
page  21.  .        ■  -  • 

A  system  of  protecting  the  temperature  of  the  column  by  minimizing 
temperature  effect  of  the  surroundings  is  a  guard  system.    The  temperature 


A 


■  (2M 
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of  the  column  is  guarded  by  the  thermostat  air  bath.    The  g\aard  system 
can  be  illustrated  by  its  electrical  analog  shown  in  Fig.  2.    The  circuit 
can  be  characterized  by  a  time  constant  T  =  R^C(l+A).    When  this  time 
constant  is  large,  compared  to  the  time  required  for  one  experiment,  the 
system  may  be  treated  as  if  the  column  were  mo-unted  adiabatically.    It  , 
was  found  experimentally  that  the  system  has  a  time  constant  of  200 
minutes,  or  approximately  20  times  the  duration  of  an  experiment. 

Although  the  basic  circuit  of  the  SCR  controller  responds  to  resis- 
tance changes,  the  current  through  the  recommended  resistance  sensor 
activates  a  safety  circuit  to  guard  against  open  and  short  circuits  at 
the  input.    It  was  found  that  the  circuit  shown  would  operate  satisfac- 
torily with  a  minimum  of  modification  of  the  SCR  \init.    The  variable 
output  current  of  the  Philbrick  amplifier  through  fixed  resistance  R 
is  equivalent  to  a  varying  resistance  in  the  SCR  circuit  and  safety 
circuit . 

The  frequency  response  of  the  preamplifier  is  limited  by  the  0.3 
I-IF  capacitor  in  the  feedback  path,  to  attenuate  any  60  cycle  noise  and 
to  eliminate  any  spurious  oscillation  of  the  overall  loop.  Maximum 
loop  gain  was  sufficient  to  drive  the  SCR  output  from  nought  to  satura- 
tion with  a  temperature  change  of  approximately  l.^C 

Heating  rate  of  the  thermostat  at  saturation  power  was  found  to  be 
greater  than  30° C  per  minute. 
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FIGURE  2 


Electrical  Analog  and  Equivalent  Circuit  of  Heat 
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DESIGN  OF  THERMOSTAT: 
II.  Construction 

In  isothermal  gas  chromatography,  the  column  should  have  a  constant 
and  uniform  temperature.    In  programmed  temperature  gas  chromatography, 
one  should  have  a  known  temperature -time  profile  on  the  column  for  a 
given  point  on  the  column  and  for  any  time  t  all  points  on  the  column 
should  be 'at  the  same  temperature.    The  thermostat  was  designed  so  that 
the  column  would  have  the  desired  temperature  behavior. 

Heat  is  generated  directly  in  the  column  by  means  of  electrical 
current,  and  the  air  bath  around  the  column  is  maintained  at  the  same 
temperature  as  the  column  by  means  of  an  auxiliary  heater.    Figure  3  is 
a  schematic  diagram  of  the  thermostat.    Air  is  circulated  downward 
around  column  container  5  through  space  1  and  upward  through  annular 
space  2,  containing  the  heater  coils  3.    The  gas  chromatographic  column  k 
is  mounted  inside  column  container  5  which  is  suspended  by  means  of  thin 
aluminum  lead-in  tubes  6.    The  lead-in  tubes  serve  both  as  containers  for 
all  electrical  wiring  connected  to  the  column  and  as  support  for  the 
column  container.    The  temperature  of  the  thermostat  is  controlled  by  a 
Honeywell  proportional  temperature  controller.  Model  R7187A.    The  con- 
troller is  driven  by  the  output  of  a  differential  thermocouple  (8)  moni- 
toring the  difference  in  the  temperature  between  the  column  and  the  air 
bath  in  the  thermostat.    The  signal  from  the  thermocouple  is  amplified 
by  a  Philbrick  SP656  amplifier. 

The  injection  port  and  heater  are  the  same  as  used  in  an  Aerograph 
(10)  gas  chromatograph;  they  are  mounted  on  the  lid  of  the  thermostat. 
The  injection  port  temperature  control  is  a  variable  transformer;  thus 
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TABIE  1  , 
Components  of  Fie,"  ire  3 

1.  Air  circulation  space 

2.  Air  circulation  space 

3.  Heaters 

k.  Gas  chromatographic  column 

5.  Column  container 

6.  Lead-in  tubes 

7.  Lid  with  insulation 

8.  Thermocouples 

9.  Sq.uirrel  cage 
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it  will  hold  constant  temperature  only  at  constant  heat  loss.    The  ideal 
injection  port  would  hold  any  temperature  desired  and  not  lose  any  heat 
to  the  column.    However,  the  injection  port  is  thermally  isolated  from 
the  column  by  mounting  it  on  the  top  of  the  thermostat.    The  connecting 
tube  between  the  column  and  the  injection  port  passes  through  the  air  bath 
where  any  needed  cooling  or  heating  occurs.    The  moving  phase,  nitrogen, 
is  preheated  in  the  air  bath  before  it  goes  to  the  injection  port. 

The  chromatographic  instrument  has  several  components  as  shown  in 
Fig.  h.    The  electrometer  of  the  Aerograph  60  A  and  its  hydrogen  ioniza- 
tion detector  are  the  detector  and  amplifier  for  the  system.     The  thermo- 
stat has  been  described.    The  read-out  system  is  a  Texas  Instruments 
recorder  (ll).    The  potentiometer  used  to  measure  temperatures  is  made  • 
by  Rubicon  and  has  a  resolution  of  0.001  millivolt.    The  readings  were 
made  to  the  nearest  0.001  millivolt. 

Figure  5  shows  the  gas  flow  diagram  for  the  chromatography  set  up. 
The  details  for  the  control  loop  are  shown  in  Fig.  6. 
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TABLE  2  ■  , 

Components  of  Figure  k 

A.  Variable  transformer 

B.  Electrometer,  (Aerograph  6o  A) 

C.  Step  down  transformer^  100  amp,  10  volt  output 

D.  Out  put  recorder  (Texas  Instruments  Servo  Recorder) 

E.  Line  voltage  120  VAC 

F.  Igniter  •  -  • 

G.  Signal  cable 

H.  Thermostat 

I.  Thermocouple  wires  (temperature  indicating) 
J.  Cold  junction  0°C 

K.  Potentiometer  • 

L.  Differential  thermocouple  wires 

M.  Power  line  -to  thermostat  heaters 

0.  Preamplifier 

P.  Honeywell  controller 

■  -.  *  . 
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Equipment  Diagram 
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FIGURE  5 


Diagram  of  Instrument  Gas  Flow 
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FIGURE  6 

Philbrick  Amplifier  and  Honeyvell  Controller 
Wiring  Diagram 
C 


I  1 

^  Fine  set  point 

A.  Coarse  set  point  with  reference  ; 

B.  Philbrick  amplifier  #  SP656 

C.  Capacitor,  0.3)aF 

D.  Power  supply,  ±  15  VDC,  6.3  VAC 

H.  Honeywell  proportional  temperature  controller  R7187A 

R-^  IM  ohm 

R^  100  ohm  '■  '  •  ■  ■ 

R^  1000  ohm 

R2j_  IM  ohm 
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OPEEATION  OF  THE  PROGRAMMED  TEMPERATURE 
GAS  CHROMATOGRAPH 

For  satisfactory  operation  of  the  programmed  gas  chromatography  the 
starting  conditions  should  be  reproducible.    The  injector  port  should  be 
hot  enough  to  q.uickly  vaporize  the  highest  boiling  fraction.    The  injec- 
tion port  must  be  at  the  same  temperature  at  the  start  of  each  run.  The 
temperature  of  the  column  must  also  be  at  the  same  starting  temperature 
if  identical  runs  are  desired.  '  If  the  starting  temperature  is  higher, 
the  peaks  will  elute  earlier  but  at  a  predictable  temperature .  The 
starting  conditions  for  the  experiments  were  kept  constant  by  allowing  the 
column  to  cool  to  a  constant  temperature.    This  temperature  was  repro- 
ducible between  runs.    The  injection  port  loses  the  heat  gained  during  the 
previous  run  and  its  temperature  is  easily  reproduced  by  allowing  con- 
ditions to  come  to  an  equilibrium.     The  amount  of  heat  the  injection 
port  gains  during  a  run  does  not  affect  the  temperature  of  the  column, 
because  there  is  no  heat  transferred  to  the  column  from  the  injection 
port.    If  the  injection  port  is  below  the  boiling  point  of  any  of  the 
components  injected,  the  heat  gained  by  the  injection  port  affects  the 
rate  of  evaporation  of  high  boiling  fractions  and  thus  their  retention 
time.    It  was  found  that  the  temperature  of  the  port  at  the  end  of  the 
run  was  reproducible  if  the  starting  conditions  were  the  same. 

For  each  run  gas  inlet  pressure  was  adjusted  so  that  V^/p^  'was  I.77. 
The  starting  temperature  of  both  the  column  and  the  injection  port  was 
checked  and  adjusted  if  needed.    The  recorder  was  started  and  when  'the 
pen  crossed  a  marker,  the  sample  was  injected  and  a  stop  watch  was 
started.    The  Honeywell  temperature  set  and  its  proportional  band  should 
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be  pre-set  at  optimum  conditions.    When  ten  seconds  have  elapsed,  power 

is  supplied  to  the  ^s  chromatographic  column.     (See    Fig.  l8  for  a  typical 

heating  curve . ) 

The  procedure  followed  for  constant  temperature  chromatography  was 
the  same  as  that  for  programmed  gas -liquid  chromatography  except  that 
ohmic  heating  of  zhe  column  was  omitted.     To  obtain  a  constant  temperature 
the  column  thermocouple  was  replaced  by  a  reference  signal  from  a  battery 
and  voltage  divider.    Thus  this  constant  signal  forced  the  guard  ring  to 
keep  an  air  bath  of  constant  temperature  around  the  column.     The  tempera- 
ture of  the  air  bath  and  the  column  is    directly  proportional  to  the  sig- 
nal from  the  battery  and  voltage  divider.     (See  Fig.  7  for  circuit  diagram.) 

I 
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PERFORMANCE 

The  linearity  of  the  heating  rate  of  the  column  and  the  predicta- 
bility of  the  retention  times  for  compounds  is  a  measure  of  the  performance 
of  the  instrument.    The  heating  characteristics  of  the  column  and  the 
linearity  of  heating  rate  were  determined.    Retention  times  and  retention 
temperatures  of  compounds  were  predicted  from  isothermal  data  and  compared 
to  the  actual  values  obtained  from  programmed  gas-liq[uid  chromatography. 

I.    Thermal  Tests 

Direct  heating  is  useful  only  if  the  temperature  distribution  along 
the  column  during  the  heating  cycle  is  unifonn.    Several  experiments  were 
carried  out  to  determine  this.    The  purpose  of  the  experiments  was  to  find 
a  means  to  sense  the  temperature  on  a  column  during  a  temperature  change 
and  to  determine  the  temperature  profile  along  a  typical  column. 

Several  methods  of  determining  the  column  temperature  were  tried. 
Thermocouples  were  the  sensors  for  all  experiments;  the  method  of  attach- 
ment to  the  column  and  the  size  of  the  thermocouple  wire  was  varied.  It 
was  desirable  to  have  a  thermocouple  which  could  quickly  respond  to  a 
small  change  in  temperature.    The  material  finally  chosen  was  wire, 
0.013  centimeter  diameter  copper  and  constantan  (B  &  S  #36).    The  time 
lag  for  these  thermocouples  in  measuring  a  cliange  in  still  air  from  800°F 
to.  100° F  WEis  5  seconds. 

To  test  the  responses  of  the  thermocouples,  various  methods  of  attach- 
ing the  thermocouples  were  compared.  '  The  wires  were  attached  to  three 
feet  of  copper  tubing  which  was  silver  soldered  to  a  brass  block  which 
served  as  a  heat  sink.    The  brass  block  was  heated,  and  heating  and 
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cooling  curves  of  the  tube  at  the  points  of  connection  were  determined. 

In  the  experiment  four  thermocouples  were  placed  around  the  tube 
equidistant  from  the  heat  sink.    One  thermocouple  was  soft  soldered  onto 
the  tube_,  two  were  attached  with  Omega  Thermcoat,  and  one  was  taped  onto 
the  tube  with  glass  tape.    The    thermocouple  that  was  soft  soldered  onto 
the  tube  had  the  smallest  temperature  lag^  thus  it  also  had  the  smallest 
time  lag.     (See  Fig.  8  for  results.)    Assuming  that  all  points  of  attach- 
ment were  at  the  same  temperature,  the  most  desirable  method  of  attach- 
ment is  soft  soldering.     This  result  is  expected  because  a  good  electri- 
cal conductor  such  as  solder  is  also  a  good  heat  conductor. 

Soft  soldering  a  thermocouple  onto  a  column  limits  the  colmn  tem- 
perature to  the  melting  point  of  solder.    Another  disadvantage  of  using 
soldered  or  brazed  junctions  is  that  additional  unwanted  heat  capacity- 
is  added  to  the  column. 

Attaching  the  thermocouples  to  the  column  by  means  of  a  spot  weld 
offered  a  logical  alternative.    The  heat  transfer  should  be  excellent 
and  a  minimum  of  additional  heat  capacity  would  be  introduced. 

An  experiment  was  done  to  compare  the  efficiency  of  the  thermo- 
couples attached  by  a  weld  to  that  of  a  thermocouple  attached  by  lead 
solder.    There  were  two  welded  thermocouples,  one  had  the  copper  wire 
1  millimeter  from  the  constantan  wire,  the  other  had  3  millimeters  dif- 
ference between  the  two  wires.     One  millimeter  was  as  close  as  the  wires 
could  be  conveniently  and  reproducibly  attached  with  a  spot  welder.  The 
third  thermocouple  was  soldered  onto  the  column.     The  soldered  thermo- 
couple was  attached  first;  then  the  other  thermocouples  were  spot 
welded  in  close  proximity  to  the  one  soldered.    The  results  of  the  ex- 
periment are  shown  in  Fig.  9.    The  temperatures  indicated  by  each  method 
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FIGURE  9 
Group  Thermocouple  Test 


TLTTie  in  Mt'nutes 
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showed  no  significant  difference  graphically.  Thus  spot  welded  thermo- 
couples will  sense  the  column  temperature,  and  3  millimeters  difference 
between  the  thermocouple  ends  introduced  negligible  error. 

The  temperature  along  a  column  must  be  uniform.     One  problem  for  a 
column  which  is  heated  ohmically  is  the  occurrence  of  hot  spots  due  to 
local  work  hardening  of  the  metal  which  increases  the  electrical  resis- 
tance.   Work-hardened  areas  due  to  bends  in  the  tubing  will  cause  hot 
spots.    An  experiment  was  done  to  find  out  if  these  effects  were  sig- 
nificant.   Two  thermocouples  were  placed  on  a  stainless  steel  l/8  inch 
O.D.  tube,  one  on  a  straight  portion  and  one  on  a  sharp  curve.    The  tube 
was  wrapped  in  glass  wool  and  heated  ohmically.    The  results  are  shown  in 
Fig.  10.    The    curved  part  of  the  tube  definitely  heated  faster  than  the 
Straight  part.    The  cooling  curve (Fig.  11^  indicated  that  both  portions 
of  the  curve  were  equally  insulated.    The  apparent  discrepancy  of  the 
temperature  between  the  two  curves  is  due  to  the  fact  that  the  two  points 
sensed  started  •cooling  at  the  same  time,  but  at  different  temperatures. 
Thus  we  can  detect  hot  spots  that  result  from  work-hardening. 

The  column  which  was  used  in  the  chromatographic  experiments  was 
likewise  tested  for  uniformity.    The  sensors  were  spot  welded  on  the 
column  and  connected  to  the  lead-out  wires.    The  copper  side  of  the  36 
gage  thennocouple  was  connected  to  a  22  gage  copper  wire.    All  wires 
are  insulated  from  one  another  but  not  from  the  column.     The  potential 
was  measured  on  a  floating  potentiometer.    The  column"  remained  isolated 
electrically.    Temperature  profile  of  the  carbowax  column  was  determined 
during  a  programmed  temperature  run,  10°c/minute.    A  total  of  ten 
thermocouples  was  placed  on  the  column,  particularly  at  suspected  hot' 
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FIGURE  10 

Test  for  Hot  Spots  Due  to  Work  Hardening: 
Heating 


8.0 


7.0  ■ 


6.0 


S.o 


o 


JO- 


2.0 


1.0- 


O  on  cuKue  Cu-'eU) 

A  TC  on  sf-t-aishtCu^U) 
□      TC  on  S+raigh-i  CUT) 

CrC  i.s  Thertr\oce>uf\e) 
(HT  is  Har</  Thermocoat) 


4  6-  -8 

Time   in  Mtnuies 


10 


31 


FIGURE  11 

Test  for  Hot  Spots  Due  to  Work  Hardening: 
Cooling 


□  TC  on  straight  (weld  and  HT) 
O  TC  on  curve  (weld) 
(HT  is  hard  Thermocoat) 
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or  cold  spots.    See  Fig.  12  for  an  illustration  of  the  placement  of  the 

thermocouples.    The  results  of  the  experiment  are  shown  in  Fig.  I3.  The 

time -temperature  plot  of  the  run  by  least  squares  method  indicated  that 

the  column  had  a  uniform  temperature  at  any  time  T.    A  simple  regression 

analysis  of  the  data  verified  the  linear  relationship  of  the  data.  The 

deviation  about  the  regression  was  O.OO165  millivolts  which  is  equivalent 

to  0.2°C.    The  variance  of  millivolts  was  0.00002  mv^  or  0.0004°C^.  The 

2  2 

variance  of  the  slope  was  0.00093(mv/min. )    or  O.OI8  (°C)  .     The  correla- 
tion coefficient  was  0. 99952.    All  the  thermocouples  indicated  the  s'ame 
temperature  at  different  points  along  the  curve  when  the  readings  were 
corrected  for  the  heating  rate.    Thus  any  end  effect  is  within  the  limits 
of  error. 

The  thermostat  was  developed  to  have  a  large  heat  transfer  resistance 
or  a  small  heat  transfer  coefficient  (explained  in  the  design  theory 
section).    The  time  constant,  T,  for  the  thermostat  is  a  measure  of  how 
well  the  column"  has  been  insulated  by  the  guard  system. 

The  time  constant  was  determined  experimentally.    The  cooling  curve 
of  the  column  was  determined.    The  Honeywell  controller  temperature  off- 
set and  the  preamplifier  offset  were  the  same  as  that  used  for  the  pro- 
grammed runs.    The  column  and  thermostat  were  heated  to  167° C  and  then 
the  system  was  allowed  to  float.    The  column  slowly  cooled  (Fig.  l4). 
The  time  constant,  T,  was  determined  to  be  the  time  for  67/0  of  the  tem- 
perature difference  to  diminish.    A  time  constant  of  200  minutes  was 
found  for  the  thermostat.  . 


FIGURE  12 
Placement  of  Sensing  Thermocouples 
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Position  indication  for  thermocouples  on  column  if  straightened 
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II .    Programmed  Temperature  Gas 
Chromatography  Tests  ; 

The  operation  of  the  thermostat  vas  tested  chromatogeraphically  by 
a  comparison  of  the  experimentally  determined  programmed  temperature 
retention  data  with  data  predicted  from  isothermal  runs  under  the  assump- 
tion of  linear  heating  rate.    Correlation  of  programmed  temperature  data 
with  isothermal  data  for  known  time -temperature  profiles  has  been  well 
established  by  others  (12). 

The  isothermal  chromatography  provides  data  to  be  used  for  pre- 
dictions of  retention  times  for  the  programmed  runs.    The  operating  con- 
ditions must  be  compatible  to  both  isothermal  and  programmed  run  pre- 
requisites.   The  isothermal  data  should  be  accurately  known  and  repro- 
ducible so  that  the  data  introduce  a  minimum  error  into  the  predictions 
of  the  retention  times  for  programmed  runs. 

Isothermal  chromatography  operating  conditions  are  constant  mass 
flow  rate  and  constant  average  velocity.    This  is  achieved  by  having 
constant  inlet  and  outlet  pressures  of  the  carrier  gas..    Most  detectors 
function  well  with  these  conditions.    However^  in  programmed  gas-liquid 
chromatography  the  mass  flow  rate  of  the  carrier  gas  changes  with  column 
temperature.    Thus  detectors  such  as  the  thermal  conductivity  detectors 
which  are  commonly  used  with  programmed  chromatography  cannot  be  used 
unless  measures  are  taken  to  maintain'  constant  mass  'flow.     Constant  mass 
flow  is  achieved  by  varying  the  inlet  and  outlet  pressures^  but  then  the 
average  gas  velocity  and  pressure  drop  are  complicated  functions  of  tem- 
perature.   Because  it  is  impossible  to  have  constant  volume  flow  in  the 
column,  the  most  convenient  operation  is  at  constant  inlet  and  outlet 
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pressures.    Under  these  conditions  the  velocity  profile  along  the  column 
is  constant,  and  only  the  average  velocity  changes  with  temperature.  A 
detector  insensitive  to  mass  flow,  the  hydrogen  flame  ion  detector,  was 
used  to  simplify  xhe  technique. 

For  isothermal  work  the  control  circuit  was  altered  to  function  with 
a  voltage  divider  and  battery.  (See  Fig.  7  for  the  circuit  details.)  No 
thermal  drift  could  be  detected  during  isothermal  runs  with  the  equipment 
used. 

The  retention  time,  t  ,  depends  upon  the  compound,  inlet  and  outlet 
R 

pressures,  temperature  of  the  column,  the  partition  coefficient,  the  amount 

of  liquid  and  gas  phase  contained  per  unit  of  the  column  and  the  length  of 

the  column.     (Refer  to  theory  section  and  equation  number  11.)  Operation 

at  constant  inlet  and  outlet  pressiire  gives  the  simplest  mathematical 

equations  for  treating  the  average  carrier  gas  velocity,  a  constant  for 

isothermal  work  and  a  function  of  temperature  only  for  programmed  gas- 

liquid  chromato'graphy.    In  both  isothermal  and  temperature  programmed 

gas  chromatography,  the  partition  coefficent  will  be  the  same  at  a  given 

temperature,  T,  which  is  known  as  a  function  of  time.    In  isothermal 

chromatography  the  velocity  at  position  x  depends  only  on  the  inlet  and 

outlet  pressures,  but  in  programmed  temperature  gas  chromatography  the 

velocity  also  depends  on  the  temperature.     The  other  paramaters  which 

affect  the  retention  time  are  the  physical  characteristics  of  the 

column;  these  do  not  change  with  time  or  temperature  and  are  constant 

for  all  runs.  ,  . 

The  isothermal  retention  time,  t  ,  at  a  specific  temperature  can  be 

r 

related  to  peak  velocity  in  linearly  programmed  temperature  experiments 
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at  that  same  temperature  if  the  time  and  temperature  can  be  accurately- 
related.  The  retention  times  for  a  series  of  compounds  at  20  tempera- 
tures were  obtained. 

The  mixture  used  for  these  experiments  consisted  of  eight  linear 
hydrocarbons:    n-hexane^  n-heptane,  n-octane,  n-nonane,  n-decane,  and 
n-dodecane.     The  percent  composition  varied  so  that  hexane  was  least  and 
dodecane  was  most.    Samples  of  O.5  1-^1  were  used  for  all  runs;  they  were 
injected  with  a  10  m  Hamilton  syringe.    Time  was  noted  on  the  strip 
chart  every  three  minutes  during  the  first  10  minutes  and  then  every 
15  minutes  for  the  rest  of  the  experiment.    When  it  was  convenient^  the 
time  of  peaks  maxima  were  noted  on  the  strip  chart.    A  typical  isothermal 
chromatogram  is  shown  in  Fig.  15- 

The  isothermal  runs  were  made  at  approximately  10° C  intervals  be- 
tween 62° C  and  200° C    The  retention  times  of  the  components  were  cal- 
culated from  the  distance  on  the  strip  chart  between  injection  and  elu- 
tion  and  detection  of  the  compounds.    This  distance  is  directly  related 
to  time,  the  exact  relation  depending  upon  chart  speed.     The  distances 
were  measured  to  the  nearest  O.O8  centimeter ,  or  0.12  seconds.    The  iso- 
thermal data  are  enumerated  in  Table  3- 

The  agreement  between  the  retention  time  calculated  from  the  chart . 
and  the  retention  time  noted  on  peaks  was  within  1^  for  all  noted  com- 
pounds with  a  retention  time  over  12  seconds.    An  analytical  relation 
between  net  retention  time  t'  and  temperature  was  not 'found.    Golay  (12) 
related  the  retention  time  to  temperature  by:  • 

■^R  =  (T^/T)Ye-A/(T-To)R  ^  (25)^ 


ho 


TABLE  3 

Compilation  of  Retention  Data 


Temp 
°A 

Net 

Retention  Time 

in  seconds 

n-c5 

n-C8 

n-Cq 

n-c^p 

n-c^^ 

.  -I.T-  

35^^.8 

9.6 

20.2 

41.5 

80.9 

179-4 

590.0 

1089.3 

2069.0 

361.7 

9-0 

18.3 

35.7 

65.9 

145.0 

433.0 

817.4 

1506.5 

369.^ 

7-1 

17.2 

32.7 

32.7 

125.6 

355.2 

640.8 

1161.2 

37^^.3 

9-6 

16.5 

29.0 

51.5 

111.9 

286.5 

461.0 

854.-7 

381.4 

9.6 

15.9 

27.1 

44.6 

97-5 

299.0 

399.0 

699.0 

390.5 

9.6 

15.0 

24.2 

4o.i 

87.5 

187-7 

319-0 

699.0 

396.2 

8.0 

13.3 

20.0 

34.0 

78.1 

155.2 

255.5 

424.3 

405.0 

7.0 

10.0 

l4.0 

21.0 

53.2 

78.3 

122.1 

191.8 

406.9 

7.0 

11.2 

18.7 

28.7 

64.4 

139-0 

187.7 

303.1 

413-5 

6.2 

9-0 

l4.6 

24.6 

58.1 

92.1 

143.7 

228.7 

424.6 

7.0 

9-0 

l4.o 

20.8 

53.1 

76.5 

119-0 

182.7 

433-0 

6.5 

9-6 

13.3 

19.6 

49.4 

63-1 

96.2 

145.8 

433-0 

6-5 

9-0 

l4.o 

18.7 

50.3 

84.7 

■  96.5 

145.8 

443.9 

6.5 

9.6 

12.1 

18.0 

45.6 

51.2 

75.5 

110.5 

455.7 

6.5 

6.5 

12.4 

17-1 

43-5 

45.0 

64.6 

92.1 

463.1 

43-7 

40.8 

54.6 

77.1 

477.0 

■ 

4o.o 

32.1 

43.3 

57.7 

477.0 

4o.o 

•  44.0 

59-0 

335.2 

25.8 

74.6 

159 

349.4 

1424.0 

2928.7 

346.0 

1 

9.6 

24.6 

92.7 

108 

233-8 

850.2 

1679.9 

3211.7 

hi 


where  t    is  retention  time,  T    is  starting  temperature  in  °K,  T  is  the 
independent  variable  in  °K,  A  and  y  are  constants.    The  constant  A  is 
the  free  energy  change;  y  is  a  viscosity  term  and  R  is  the  molar  gas  con- 
stant.   Equation  (25)  did  not  fit  the  curve  within  1^  error.    To  relate 
net  retention  time,  T'^,  and  temperature,  a  graph  was  made.    A  typical 
example  is  Fig.  l6.    The  net  retention  time  was  related  to  temperature 
by  the  curve.    The  graph  was  used  to  obtain  the  net  retention  time  at 
intervals  of  1°  from  335°K.    These  data  were  used  to  predict  the  reten- 
tion time  of  the  compound  at  a  specific  linear  time -temperature  program. 

The  programmed  temperature  conditions  have  already  been  stated;  they 
were  selected  to  be  compatible  with  the  isothermal  runs.    The  method  of 
injection  was  the  same  as  that  used  in  isothermal  work,  except  that  the 
power  to  the  column  was  switched  on  10  seconds  after  injection.  The 
chart  speed  was  I.5  inches  per  minute  throughout  the  experiment.  The 
distances  between  injection  and  detection  were  related  to  retention  time. 
A  typical  programmed  temperature  chromatogram  is  shown  in  Fig.  I7. 

The  temperature  of  the  column  during  the  program  was  measured  as  a 
function  of  time.    From  the  data  the  heating  rate  of  column  was  proven  to 
be  linear  by  linear  regression  analysis.    A  typical  heating  curve  is  shown 
in  Fig.  18.    For  carefully  prepared  programmed  heating  rates,  the  sum 
about  the  regression  was  0.177°C,  the, standard  deviation  of  the  intercept 
was  0.2°c/min.     The  data  for  the  programmed  temperature  runs  will  be 
given  later  in  this  section. 

The  column  temperature  should  rise  immediately  when  ohmic  heating  is 
started.  The  ideal  thermostat  will  have  an  instant  change  to  the  heating 
rate  desired  when  the  program  is  started.  "  However,  the  temperature*  of  the 
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FIGURE  18 


Heating  Curve  for  Ideal  Time -Temperature  Program 


1^5 

column  did  not  change  instantly.    The  thermostat  and  column  used  could  be 
operated  so  that  there  was  a  temperature  lag  or  lead,  neither  of  which  is 
desirable  for  linear  programs.    The  thermostat  controls  were  set  for  a 
temperature  lead  and  then  altered  so  that  the  column  would  border  on  the 
conditions  that  would  allow  a  temperature  lag.    Representative  results  are 
shown  in  Hahle  k.    The  column  could  be  operated  so  that  the  time  differ- 
ence between  actual  start  and  calculated  start  was  as  short  as  O.5 
seconds  or  less. 

Computation  of  programmed  temperature  retention  times  from  isothermal 
data  were  based  on  the  relationship 

r^R 

L  =  I      u  rtt  (26) 
s 

in  which  u^  was  the  instantaneous  sample  velocity.    Three  factors  entered 
into  the  calculation  of  u^:    position,  average  gas  flow  rate,  and  parti- 
tion coefficient.    The  position  factor  was  computed  from  equation  (5),  the 
other  two  factors  were  obtained  from  isothermal  data.    An  I.B.M.  l6kO  com- 
puter was  programmed  to  compute  the  progression  -of  a  sample  component  in 
a  series  of  small  time  increments,  until  the  total  linear  distance 
traveled  exceeded  column  length.    The  total  time  required  was  determined 
from  an  average  of  time  required  for  a  component  to  reach  the  approximate 
end  of  the  column  and  the  next  time  period  during  which  the  component 
e lutes  from  the  column. 

■  In  order  to  solve  equation  (26)  4:he  right  side  of  the  equation  must 
be  a  function  of  t  only.    The  instantaneous  sample  velocity  is  a  function 
of  k  and  u^  and  position 


TABLE  k 
Representative  Results 


Plateau 
time  delay- 

Heating 
rate 

Difference  between  actual  ' 
delay  and  calculated  delay 

sec. 

°  c/min 

sec. 

10 

k.90 

-k.k* 

-O.h 

10 

9-19 

-2.5* 

-o.h 

10 

13.78 

-O.k 

-0.1 

10 

5.02 

0.0 

0.0 

20 

12.76 

0.5 

0.1 

*  The  large  delay  vas  due  to  operator  error. 


Therefore 


s         ■   '  + 

^r  _  isothermal 

Thus  equation  (26)  can  be  written 


u    =  -^-^  •  27) 

^  1+k 


u    =  f(x)-^   •  (28) 


L  =  f(x)^(dt).  ^  (29) 


Where 


r  isothermal 


^r  isothermal      ■^(^)'  (30) 


Now 


^  L 


L  =  E  f(x);        (At  )=  E  f(x)f(t)L  At^.  (31) 
.  i=0         riso.  i=o 

Where  nAt  =  t^,  and  f(x)  can  he  computed  from  constants       and  p^. 

Equation  (9)  relates  average  retention  time  to  velocity  at  the  outlet 

for  the  moving  phase.    Because  u^/u  can  be  substituted  for  p/p^^  equation 

(8)  may  be  changed  to      ■  " 

P    r  (P./P  )   -1  -1 
,1'  o' 

where  u^  is  velocity  at  the  end  of  the  colLimh;  u  is  the  average  velocity^ 
p^  and  p^  are  inlet  and  outlet  pressures  respectively.    Equation  (32)  can 
be  simplified  because  V^/v^  always  equals  1.77  and  the  column  length  is 
constant :  -  ■  ■  ■  ' 

■     u   =  ^14—   •  •  (33) 

°     t^  /1.4207  ^^^^ 


k8 


where  t'^^  is  net  retention  temperature  for  a  compound  at  temperature  T. 
Thus  it  is  convenient  to  define  a  t''  as: 

I 

Then  t''^^  is  a  function  of  net  retention  time  and  can  be  used  to  relate  u 
to  the  velocity  of  the  peaks  at  the  end  of  the  column^  u^. 

The  equation  used  for  prediction  of  the  velocity  of  a  peak  at  position 
X  at  time  t '  is : 

-  /  (35) 

yf[(p./p,)2_i)+(p./p^)2  . 

Equation  35  can  be  rewritten  as  : 

61 


u  = 

X 


*R  7-0.03496x^-3. 1329 


(36) 


Equation  (36)  shows  how  t'     can  be  related  to  u  .    This  equation  can 

R  X 

be  substituted  into  equation  (26)  and  then  equation  (26)  can  be  integrated. 
This  problem  was  done  on  the  IBM  l4lO  computer  by  an  average  step  inte- 
gration.    The  form  of  the  computer  program  is: 

61 


J    f (u^)dt    E  — — :       ■   .  At'  (37) 

•Jq  i=o      V-0.03496X.+3.1329  ^""^ 


Where  t '    -  nAt . 
R 


The  computer  equation  was: 

61 


=  (38) 


where  x  is  zero  when  the  component  was  injected.    The  t"^  was  a  function 
of  both  the  component  and  the  temperature  and  was  obtained  from  isothermal 
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data.    The  DI^^  vhich  equals  X;,  the  position  of  the  peak  (see  Fig.  19); 
was  zero  when  x  is  zero  and  computable  when  x  was  greater  than  zero. ■ 
In  equation  (38) 


p  =  K(Pi/pQ)^-i  =  0.03496] 


(39) 


and 


Q  =  (P^/Pq)^  =  3.1329  .  (40) 

When  X  was  11^  the  time  elapsed  was  calculated  by  iR+-10  (R  was  inverse 
heating  rate  in  seconds/^ C). 

The  best  velocity  for  a  peak  during  an  increment  At^  is  the  average 
velocity  of  the  gas  obtained  due  to  the  change  in  position  and  tempera- 
ture.   This  average  was  approximated  by  using  the  u^  to  predict  where  the 
gas  will  be  at  (x+l)  RflO^  then  by  finding  the  velocity  there  by 

61 


~  ""'^c^i  V-P(Dl^i)+Q 


(^1) 


'■xt-l^ 

where 


At(u  )  =  ADI'  ^  .  (42) 

^  X 

The  ADI'^  was  the  predicted  distance  the  peak  would  travel  during  the 
time  increment.  Then 

\  =  -^T^ '    .  ■     ..  ^^^^ 

The  step  integration  was  done  by  finding  the  distance  traversed  during 
time,  At,  by: 
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FIGURE  19 


Diagram  of  Method  for  Calculation  of 


fc 


DI. 


DI 


fo 


I  -  =  ^■f"^'^^^'*^^  stop    wken     su/*M'Matione««^f5  L, 

~  J^S''  ^  z  nit  tiO 

o  Is 
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where  t  is  equal  to  R.     The  distance  traveled  during  time  xRMO  was 

DI^=  Z    ADI-.  (i|6) 
i=l  i 

The  net  retention  time  for  a  compound  at  a  given  temperature  represented 
one  piece  of  information.      The.  pieces  of  information  were  arranged  into 
a  progression  so  that  the  retention  time,  at  1°C  increase,  followed  the 
previous  one.    This  information  was  used  in  the  computation  and  referred 
to  as  the  "data  pack". 

The  net  retention  times  from  programmed  runs  were  compared  with 
the  calculated  values  (Table  5)-    The  time  required  for  the  temperature 
to  rise  1°C  after  the  plateau  was  R.    The  rate  of  temperature- rise,  R,  was 
the  value  used  in  calculating  the  net  retention  time.     The  deviation, 
abbreviated  dev.  in  the  table,  was  the  difference  between  observed 
value  and  the  calculated  value. 

The  percent  deviation  would  be  expected  to  decrease  as  the  reten- 
tion times  increase,  if  the  heating  rate  had  random  errors.     (The  sum  of 
the  deviations  would  be  zero.)    However,  the  observed  time  deviation  was 
nearly  constant  and  always  positive  so  the  error  was  not  pure  random 
heating  error. 

If  the  system  had  a  time  lag,  the  largest  percent  error  and  smallest 
time  deviation  would  be  for  the  least  retained  component.    The'  least 
retained  component,  hexane,  which  was  mobile'  at  starting  conditions,  had 
similar  time  deviations  as  the  more  retained  components  which  were  im- 
mobile at  the  starting  conditions.    Also  the  same  time  deviations  were 
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TABLE  5 


Comparison  of  Experimental  Values  With  Calculated  Values 


Component 

R 

t 

R 

Recalculated  t'-^ 

obs 

cal 

dev 

dev 

sec/°C 

sec 

sec 

sec 

io 

sec 

% 

Hexane 

12.0 

12.7 

9.5 

3.2 

34.8 

0.4 

4.2 

Heptane 

It 

25.2 

22.7 

2.5 

11.0 

-0.3 

-1-3 

Octane 

II 

50.6 

48.5 

2.2 

4.5 

-0.6 

-1.2 

Nonane 

tt 

92.7 

89.6 

3.1 

3.4 

0.3 

0.3 

Decanet 

If 

165.2 

179.8 

-l4.6t 

-8.it 





Dode  cane 

It 

398.9 

397.0 

1.9 

0.5 

-1.1 

-0.3 

Hexane 

7.02 

14.6 

11.3 

3.3 

29.2 

0.5 

4.4 

Heptane 

ri 

30.2 

26.8 

3-4 

12.7 

0.6 

2.2 

Octane 

.. 

58.9 

57.6 

1-3 

2.4 

-1.5 

-2.6 

Nonane 

11 

106.4 

102.5 

3.9 

3.8 

1.1 

1.1 

Decanet 

tt 

170.8 

187.8 

-17. ot 

-l.lt 

Dode  cane 

11 

348.9 

358.9 

-10.0 

-2.8 

-12.8 

-3.6 

Hexane* 

4.70 

12.7 

9.6 

3.1 

32.2 

0.3 

3.1 

Heptane* 

II 

23.3 

20.0 

3.2 

16.0 

0.4 

2.0 

Octane* 

If 

43.9 

4o.i 

3.8 

9.5 

1.0 

2.5 

Nonane* 

" 

77.1 

70.2 

6.9  • 

9.8 

4.1 

5-9  ' 

Decanet* 

123.9 

133.5 

-9.6t 

-7.2t 

Dode  cane* 

II 

248.9 
 .. 

242.9 

6.0 

2.5 

3.2 

1.3 

t    Deviation  for  compound  constantly  low. 


.*    The  heating  rate  is  nonlinear. 


53 


observed  for  Rl  and  R2;  this  was  not  expected.  The  time  lag  depends  on 
the  heating  rate;  therefore^  the  time  lag  should  vary  with  R.  The  error 
is  constant  for  most  compounds  in  Rl  and  R2;  hence  the  error  is  probably 
not  due  to  a  lag. 

Nonlinear  heating  would  cause  an  increasing  or  a  decreasing  devia- 
tion as  the  time  increased.    Also  the  deviations  would  be  all  positive  or 
all  negative  in  all  cases.    This  does  not  seem  to  be  the  case;  besides,  the 
deviations  in  Rl  and  R2  are  not  similar  to  the  deviations  in  R3  which  has 
a  known  nonlinear  heating  rate. 

Another  possibility  is  that  the  method  of  prediction  has  some  error. 

The  constant  time  deviation  seems  to  indicate  that  some  constant  factor 

such  as  t  ,  dead  time,  has  been  neglected.    If  the  average  deviation  is 
a 

calculated  (excluding  the  one  constantly  low  value)  and  if  this  average 
deviation  in  seconds  is  added  to  each  of  the  calculated  net  retention 
times,  the  recalculated  net  retention  time  deviations  can  be  found.  If 
the  errors  in  the  method  are  random,  both  positive  and  negative  devia- 
tions will  occur.    If  the  heating  rate  is  linear,  as  the  time  temperature 
plot  indicates,  the  percent  errors  should  be  about  equal.    The  adjusted 
net  retention  deviation  behaves  the  way  one  would  expect  and  this  indi- 
cates that  the  method  of  calculation  has  a  constant  and  removable  error. 
If  the  program  is  nonlinear,  as  is  R3;  the  corrections  will  not  help  and 
the  errors  will  not  subtract  out. 

The  average  constant  deviation  in  seconds  is  equal  to  twice  the 
error  introduced  by  the  read-out  system;  therefore,  the  least  one  can  say 
for  the  system  is  that  it  is  as  good  or  better  than  the  present  systems 
and  possibly  much  better. 


'SUMMARY  MD  CONCLUSIONS 

A  thermostat  which  has  a  predictable  linear  heating  rate  has  been 
designed  and  developed.    Solutions  to  problems  such  as  measuring  the 
column  temperatures  at  a  small  time  lag  have  been  developed.    The  major 
problems  in  ohmic  heating,  nonlinear  heating,  and  practical  insulation 
of  the  column  from  other  parts  of  the  instrument  have  been  solved  for  the 
system  used.    It  was  found  that  the  prepacked  stainless  steel  column  used 
did  not  have  any  observable  end  effects.    The  instrument  provided  useful 
data  for  both  isothermal  and  programmed  temperature  gas-liquid 
chromatography. 

This  thermostat,  which  was  made  from  readily  available  materials, 
functions  as  well  as  (and  possibly  better  than)  the  read-out  system. 
The  results  of  the  experiments  with  the  thermostat  indicates  that  a 
thermostat  for  packed  columns  and  capillary  columns  can  be  built  which 
will  have  a  predictable  linear  heating  rate  and  a  time  lag  of  3  seconds  or 
less. 

SUGGESTIONS 

The  results  of  the  experiments  indicate  improvements  that  could  be 
made  for  a  thermostat  based  on  a  guard  system.    The  injection  port  should 
be  mounted  so  that  it  is  thermally  insulated"  from  the  thermostat,  and:  the 
temperature _ effects  from  the  heat'  of  the  injection  port  and  tubes  connect- 
ing with  the  column  should  be  minimized. 

The  injection  port  should  hold  a' constant  temperature;  thus  a 'more 
elaborate  system  for  temperature  than,  a  variable  transformer  is  desirable. 
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The  time  constant  for  the  thermostat  could  be  increased  by  closer 
control  of  the  air  bath  temperature.    The  coltunn  container  will  have  to 
be  altered  to  achieve  this. 

Befoy-e  the  thermostat  can  be  operated  without  excessive  down  time, 
the  temperature  sensors  must  be  improved.    The  sensors  used  were  available, 
inexpensive,  and  functionable,  but  they  were  the  main  factor  causing 
down  time.    All  other  parts  were  satisfactorily  reliable. 
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Abstract 


A  thermostat  has  been  developed  for  use  in  programmed  temperature 
gas  chromatography.    The  design  is  based  on  direct  ohmic  heating  of  the 
column  under  conditions  of  minimum  heat  transfer  from  and  to  the  sur- 
roundings.   This  is  achieved  by  control  of  the  thermostat  air  in  such  a 
way  that  the  temperature  difference  between  the  column  and  the  thermo- 
stat air  is  kept  small  at  all  times. 

Constructions  were  developed  to  insulate  the  column  electrically 
without  introduction  of  excessive  dead  volume,  and  to  minimize  thermal 
end  effects.    Heat  effects  due  to  variations  in  electrical  resistance 
of  the  column  caused  by  work  hardening  were  foiand  to  be  negligible  if 
suitable  precautions  were  taken  during  forming  and  mounting  of  the 
stainless  steel  column.    Response  of  the  system  was  studied  by  means 
of  0.13  millimeter  copper- constantan  thermocouples  mounted  on  the 
column  by  spotwelding,  cementing,  and  soldering.    Spot  welded  sensors 
were  found  to  give  the  best  results. 

Gas  chromatographic  performance  was  studied  by  a  comparison  of 
retention  data  obtained  during  an  actual  programmed  temperature  run, 
which  was  predicted  to  be  linear  from  the  heat  input,  with  data  computed 
from  a  series  of  isothermal  runs  on  the  same  column.    The  experimental  data 
for  a  series  of  n-alkanes  agreed  to  about  2/o"  with  data  computed  under' the 
assumption  of  a  linear  time -temperature  program. 


